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Variational method has been applied to estimate Landau-Pomeranchuk-Migdal (LPM) effects on 
virtual photon emission from the quark gluon plasma as a function of photon mass. The variational 
method was well tested for the LPM effects in real photon emission. For virtual photons, LPM 
effects arising from multiple scatterings of quarks in the plasma are determined by the integral 
equations for the transverse vector function (f(pj_)) and the longitudinal function (<?(i>x)). We 
extended the variational method to solve these transverse and longitudinal equations for a variable 
set {po,9o,Q 2 }, considering bremsstrahlung and aws processes. We solved these equations, also 
by the self consistent iterations for comparing with the results of variational method. In order 
to estimate the variational parameter, we obtained empirical fits for the peak positions of the 
Px ■ f (px), Px<?(px) distributions from iteration method. We propose that the optimized variational 
parameter for virtual photon emission is approximately equal to these empirical peak position values. 
The detailed study showed that the variational method gives reliable results for LPM effects on 
virtual photon emission for photon virtuality of the order Q 2 /T 2 < 100. At low Q 2 , the peak 
positions for p± distributions of transverse vector functions for virtual photons nearly coincide with 
the peak positions of corresponding distributions of real photons. We calculated imaginary part of 
photon retarded polarization tensor as a function of Q 2 /T 2 using empirical variational parameters. 
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Study of the physical processes in quark matter as 
compared to that in hadronic matter plays a crucial role 
for identifying the Quark gluon plasma (QGP) state, ex- 
pected to be formed in the relativistic heavy ion colli- 
sions. In this context, electromagnetic processes such as 
photons and dilepton emission are important as signals 
HI 00 to identify this de-confined state. In depth study 
of photon emission processes in quark-gluon plasma were 
presented 0, . In the hard thermal loops ! 6j] (HTL) ef- 
fective theory, Compton scattering and quark-antiquark 
annihilation processes contribute at one loop level to pho- 
ton emission in quark matter. At the two loop level, the 
processes of bremsstrahlung and a crossed process of 
annihilation with scattering called aws [3, |9( contribute 
to photon emission. Importantly, these two loop pro- 
cesses contribute at the leading order 0(aa s ) resulting 
from a special effect called the collinear singularity that is 
regularized by the effective thermal masses. Owing to the 
same reasons, higher loop multiple scatterings having a 
ladder topology also contribute at the same order jlOl. fll| 
giving a decoherent correction to the two loop processes. 
These rescatterings have been resummed |lfl [ill, ef- 
fectively implementing the Landau-Pomeranchuk-Migdal 
(LPM) effects [H El Ei LPM effects arise due to 
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rescattering of quarks in the medium during photon for- 
mation time. The rescattering corrections strongly mod- 
ify the two loop contributions for bremsstrahlung and 
aws processes for real photon emission ^| . 



Photon production rates from bremsstrahlung and 
aws processes including the LPM effects are estimated 
by using Eq^ in terms of a transverse vector function 
Px ■ 5tf (p_l) E3- The resummation of multiple scatter- 
ings leads to the AMY integral equation for the trans- 
verse vector function for real photons, given in Eq|2 11]. 
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The function Kf (pj_) in Eq^ which consists of the LPM 
effects, can be taken as transverse momentum vector 
(p±) times a scalar function of transverse momentum 
pj_. The tilde sign ~ represents dimensionless quanti- 
ties in units of Debye mass mo as defined in [Tl|. The 
function p^ • 5ftf (pj_) is determined by the AMY equa- 
tion (Eq[2J in terms of energy denominator 5E given in 
Eq|3]and the collision kernel (C(£±)). We reported that 
the complex LPM effects can be very well reproduced by 
introducing the photon emission function g(x) defined 
in EqQ] of a dynamical variable x defined in Eq|S] 15]. 
In terms of a single variable function g(x), the photon 
emission rates are estimated using Eas l4l7l for any quark 
momentum, photon energy and plasma temperature. 



VARIATIONAL METHOD FOR EMISSION 
OF REAL PHOTONS 



AMY integral equation consisting of LPM effects was 
solved using the variational a ppr oach and the photon 
emission rates were reported in . In the present work, 
we follow this variational method and expand the real 
and imaginary parts of f (pj_) in terms of a basis set {<fij} 
of trial functions as shown in Eas l8llll The dimensions 
of the spaces for real and imaginary parts are N r and Ni. 
The energy function SE mn and the quantities 

Sm,T are 

calculated by the overlap integrals with the basis trial 
functions as shown in Ecis ll2ll3l The overlap integrals 
C r mn that involve the collision kernels are shown in Eall4l 
Another quantity C l mn involving imaginary part is simi- 
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Oj,T) &j,T m l ne Eas l8l9l are the expansion coefficients for 
real and imaginary parts of f(p , ) respectively. Here, the 
subscripts T and also in S m ^T represent transverse vector 
function f(pj_). Using the trial functions, some of these 
quantities can be evaluated analytically or simplified as 
shown in Eas ll5lll?l The choice of dimensions for real and 
imaginary parts should be pragmatically large enough for 
performing calculations. For real photon emission calcu- 
lations, we used two flavors, three colors, a s =0.20 and 
N r = Ni = 8. A few cases were verified for convergence 
using N r = A^=10 and 12. These results for photon 
emission rates were reported earlier |l5j . However, it is 
instructive to present the details of variational calcula- 
tions and analysis of results for the case of real photons. 
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In the variational approach, the trial functions consist of 
a variational parameter A which needs to be optimized 
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to obtain correct results. The p± • 3?f(pj_) distributions 
are usually sharply peaked and therefore the optimized 
value of the variational parameter A V (A = A 2 ,) should 
be around the peak positions of these distributions. In 
our earlier paper, we simplified this variational approach 
and extended this method to finite baryon density case 
|16| . It was shown that the variational parameter can be 
taken as A = A 2 (p ,q ,T), together with A v (p Q ,q ) — 

\po(qo + Po)/qo\°' 31 and constraint A v > 0.31. In the cal- 
culations that follow, we used this empirical result for the 
variational parameter. In order to test the empirical val- 
ues of the variational parameter, we transform the inner 
product integral p± • 9£f(p_i_) defined over [0,oo] to a fi- 
nite interval of [0,1] by defining a variable y in Eg 1201 (see 
In terms of y, this inner product is given in Eo l21l 
Then we estimate the value of p± where the integrand 
in Eq[2] peaks. If the value of the variational parameter 
A(= A 2 ) used is correct, the result should be y = | at 
this peak position p±. Therefore, finding the optimized 
value of the variational parameter is equivalent to finding 
the peak positions of these p± • 5fcf(p_i_) distributions. 
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We observed that these distributions are not very sen- 
sitive to the exact value of the A v parameter provided 
the dimension is sufficiently high. Figure. ^ shows y 
at the peak positions of these p_i_ • 5tf(p_i_) distributions. 
Figure. Of a) shows y for bremsstrahlung and IHb) for 
aws processes, as a function of \po{qo + Po)/qo\ ' ■ The 
figure includes variational calculations for a set of 2304 
values of {po,qo,T} (=24x24x4 values i.e., 24 for po, 
24 for qo and 4 temperatures). It can be seen that 
the y value is close to | for both bremsstrahlung and 
aws processes. However deviation of y from 1/2 in- 
creases with increasing temperature, as temperature de- 
pendence for variational parameter was not considered. 
The exact peak position values of these distributions 
are shown in Figure|2{a) for bremsstrahlung and Efb) 
for aws processes, as a function of \po(qo +Po)/<lo\ ' ■ 
As shown in figure, the peak positions are linear and 
the temperature dependence is not very strong. How- 
ever, when these peak positions are plotted versus a^' 31 
with x\ = |koPo(<7o + Po)/ (Tqo)\, all the 2304 points of 
{Po,Qo,T} merge for bremsstrahlung process and simi- 
larly for aws process as shown in Figs|3fa,b). Further, 
when the bremsstrahlung and aws peak positions are 
both together plotted versus x?" 31 , all the data points 
of Fig|3Ja,b) overlap as shown in Fig|^c). Therefore, 
the peak positions of these two processes follow scaling 
with x\ variable. This is not surprising because we have 
shown earlier that x\ is an exact dynamical scale for real 
photons We have fitted the data by a suitable func- 
tion and the parameters are shown in Fig|3fc). Therefore 
for real photon emission, the optimized variational pa- 
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FIG. 1: The y function constructed from the peak positions of 
real photon px • 5Rf (p±) distributions. Symbols represent the 
peak positions of pi distributions from variational method 
for different {po,qo,T} values. The variational method used 
an empirical formula for variational parameter. The y value 
should be 1/2 for truly optimized variational parameter. Fig- 
ure shows that the empirical variational parameter values 
used for real photons satisfied this constraint approximately. 



ramctcr is approximately given by the empirical fit values 
shown in Fig|3Jc) for any temperature, quark momentum 
and photon energy values. 



II. VARIATIONAL METHOD FOR VIRTUAL 
PHOTON EMISSION 

Processes that contribute to virtual photon emission 
in QGP at aa s order 01 an d the higher order correc- 
tions 0| were reported. The HTL one loop processes 
qq — > (?7*, and qg — > qj* contribute to photon polariza- 
tion tensor at the order aa s . Photon emission from QGP 
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FIG. 2: The peak positions of real photon px -5Rf(px) distri- 
butions versus (peak positions from) empirical formula. Sym- 
bols represent the peak positions from variational method for 
different {po,qo,T} values. 



as a function of photon mass considering LPM effects 
was also reported It was shown that the multiple 

scatterings modify the imaginary part of self energy as 
a function of photon energy and momentum both, espe- 
cially modifying the tree level threshold, Q 2 = 4M^ . 
The dilepton emission rates are estimated in terms of the 
imaginary part of retarded photon polarization tensor for 
virtual photons, given by 
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3IT^ M ((5) is determined by transverse and longitudinal 
functions represented by f (p_i_), g{p±) resulting from two 
integral equations. For the case of virtual photon emis- 
sion, the transverse vector function is determined by the 
Eg ESI and the energy transfer SE((p±,po, <?o, Q 2 ) 
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FIG. 3: The peak positions of real photon px • Kf(px) dis- 
tributions versus (peak positions from) empirical formula val- 
ues, for (a) bremsstrahlung and for (b) aws. Symbols repre- 
sent the peak positions of px distributions from variational 
method for different {po, qo,T} values. Figure shows that the 
x\ variable defined in text (Eg I43H for real photons is a good 
scale. The data of (a) and (b) are all plotted in figure (c). 
The empirical fit function and its parameters are shown in 
figure (c) which can be used to fix the optimized variational 
parameter. 
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For Q 2 > AMI 



For virtual photons, the coupling of quarks to longitu- 
dinal mode must be considered. This results in a scalar 
function of px which is determined by the AGMZ inte- 
gral equation given in Eo l24l with the collision kernel from 
|20j . For the case of massive photon emission, this energy 
denominator SE is modified from that of real photons by 
► M eS = Ml + ^Poro as in EqEH 
this M e g can vanish or even become 
negative. SE energy denominator can be interpreted as 
inverse formation time of the photon, which acquires de- 
pendence on photon mass in addition to the dependence 
on photon energy, quark momentum of a real photon. 

2 P± = i5E{p ± ,p ,q ,Q 2 ){{p ± ) + g 2 C F T 

1 0y2 C (^)[ f (P^)- f (P^+^)] ( 23 ) 
Po ,qo,Q 2 )g(p±)+g 2 C F T 
C(e x )[g(p x )-g(p ± +l ± )} (24) 



for the longitudinal part could be written as, 
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In above equations, r$ = po + qo and f, g are actually 
functions of Po,qo,Q represented as, f (p±,po, qo, Q 2 ), 
g{p±,Pa, qo,Q 2 ). Eg 1231 and Eq[21are identical except for 
SE energy factor. Further, Eg 1231 and the Eq[^]are sim- 
ilar on the right side of the equations, howev er the left 
side of AGMZ equation is a constant yfpo^oj- Aurenche 
et. al, solved these equations, based on a method of 
impact parameter representation |l9j . As shown in pre- 
vious section, we have solved the AMY equation for real 
photons by the variational method and a new method 
called iterations method by formulating these equations 
in terms of tilded variables. Therefore, we will trans- 
form the two above equations Eos l23l24l to tilded quanti- 
ties, and for details see [Til. l2l| . The transformation for 
px, f (p x ), SE(p ± ), C(p x ) are given by Eq[H 
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In above K e g- = We will divide the above Eg 1241 

by nip in order to get the following equation, where ab- 
sorbing 1/rriD factor, g is now re-defined. The equation 
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In the above equation, g transforms as £f(px) = 
^W-g(p±), similar to f(px) function. It implies that the 
g is larger than g by a factor rau- Therefore, when the 
solution for this Eq[2Hl is obtained and integrated over 
/ pjrp ^ff(PJ-); the result will be larger than the true re- 
sult from Eq|2]by exactly mjj factor ,21]. In the present 
work, we have solved the above EqlJHlby iterations and 
variational method with an aim to test the validity of the 
variational method for a virtual photon case. The vari- 
ational method for longitudinal part has been rc-dcrived 
as shown in Eos l3QI4l"l 22| . The basis for expansion of the 
f(i>x) is same as in Eos 181 111 The equations for SE mn 
and S m ^T remain same as in Eas ll5llo1 except for the 
change k — > k & q. 
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The equations for expansion of g(p±) and the basis func- 
tions are similar to transverse part and the equations are 
given in Eas l3QI3lfl The energy denominator and other 
quantities are shown in Eas l35l4*Tl It is important to note 
that the variational parameter A in Eas l3QI4*Tl is in gen- 
eral different for transverse and longitudinal parts. Use 
of same symbol is misleading and further, the variational 
parameters for these two parts need to be optimized in- 
dependently. This is because, in addition to overall p\ 
factor difference in • JJf(pj_), $lg(p±), the equations 
differ also in the interference terms in Eas ll9l4T1 In all 
the calculations that follow, we have used two flavors, 
three colors, a s =0.30, N r = N t = 10. and T=1.0GeV. 
Following the iterations and variational methods, we ob- 
tained the p_L distributions for the bremsstrahlung and 
aws cases for both transverse and longitudinal compo- 
nents. We obtained 350 distributions (5 for go, 10 for Q 2 
and 7 for po values) for transverse and 350 distributions 
for longitudinal parts of bremsstrahlung. Similarly the 
distributions were obtained for aws case. All these dis- 
tributions were generated in iterations method and many 
of the cases also by variational method. The results of 
these two methods were compared. The variational pa- 
rameter has been varied to optimize the distributions. It 
was observed that for a given go, the variational method 
does not give satisfactory results for increasing Q 2 values 
and further variational distributions showed large oscil- 
lations. However, when the Q 2 is low, the agreement of 
results from variational and iterations methods is very 
good. Therefore, variational method can be used for low 
Q 2 values very reliably. The iterations gave correct con- 
verging results for all go , Q 2 values studied. Therefore in 
the present work, the iterations method is taken as refer- 
ence standard for these distributions. These details 
will be presented in the next section. 



III. EMPIRICAL ANALYSIS OF THE 
SOLUTIONS OF AMY, AGMZ INTEGRAL 
EQUATIONS 

Figure 0] shows the p_i_ distributions for p± • Kf(p^) 
(transverse part of) bremsstrahlung process for five val- 
ues of qo — 50,30,20,10,5 GeV. For each g , we have 
shown distributions (in different colors in figure) for two 
different photon and quark momenta {Q 2 ,po) values. 
The symbols represent results from iterations method 
and the curves from the variational method. In many 
cases, the curves are not visible as the symbols are over- 
written on curves, suggesting that the agreement is very 
good. However, as the Q 2 increases, these distributions 



from the variational method show oscillations and devi- 
ate from the iterations method. The agreement of vari- 
ational results with iterations is good for all cases of 
Q 2 less than the higher Q 2 shown in the figures. For 
example, in the Figure ^ the deviation increases for 
any Q 2 > (50 2 - 47 2 )GeU 2 for q = 50GeV ; Q 2 > 
(30 2 -27 2 )GeU 2 forg = 30GeV ; Q 2 > (20 2 -17 2 )GeV 2 
for q = 20GeV ; Q 2 > (10 2 - 5 2 )GeU 2 for q = lOGeV 
; Q 2 > (5 2 - l 2 )GeU 2 for q = 5GeV ; 

Figure [3] shows the p^ distributions p_i_ • 5ftf (p±) of aws 
process for five values of go = 50, 30, 20, 10, 5 GeV. For 
each go, similar to previous figure, distributions for two 
quark momenta values are shown . Figure [5] shows the 
p_L distributions $tg(p±) of bremsstrahlung process for 
five values of go. Figure [7| shows the p± distributions 
for $lg(p±) of (longitudinal part of) aws process for five 
values of go- 

The p± distributions presented in Figures EE proved the 
validity of variational method. For low Q 2 this method 
can be applied to study the LPM effects in virtual photon 
emission. However, the variational parameter needs to be 
optimized as it is not known for virtual photon case. For 
the case of real photons, we have shown an empirical es- 
timate for the optimized variational parameter in terms 
of real photon dynamical variable x\. For the case of 
virtual photon emission such a simple empirical formula 
in terms of x\ is not valid and this should be a function 
also of Q 2 . Therefore, in order to predict the optimized 
variational parameter for virtual photon emission case, 
we propose to examine the peak positions values of pj^ 
distributions from iterations method. In Eqs02EHl wc 
define four dimensionless variables used in the following 
work. Especially, x\ of Eg 1451 is the relevant dynami- 
cal variable for virtual photon emission at high Q 2 and 
the variable Xi is for real photons, (inverse of x variable 
used in 0]). We searched for the peak position values 
of p^ distributions from the iterations method. These 
are functions of Po,qo,Q 2 and we searched for dynami- 
cal variables that could represent the peak positions for 
virtual photon case. 
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Figure [Hfa) shows the peak positions of the (transverse) 
p_L ■ 5Rf(p±) and (longitudinal) p±^Rg(p±) distributions 
for bremsstrahlung process at Q 2 > lOGeV 2 . The data 
is generated by solving the integral equations using it- 
erations method for the set of values {po,qo,Q 2 }- For 
the longitudinal part, the peak positions do not ex- 
ist for several values of Po,qo,Q 2 and the correspond- 
ing Kg(p^) distributions exhibit a Woods-Saxon form 
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FIG. 4: The px distributions px -JRf(px) for transverse part 
of bremsstrahlung photon emission. The symbols represent 
results of iterations. About 350 cases of {po,Qo,Q 2 } values 
were studied and some are shown in the figure and see text 
for details. Different colored symbols represent different quark 
momenta. The results of variational method are shown by 
curves. The terms It. and Var. in figure labels refer to 
iterations and variational methods. 



FIG. 5: The px distributions px • Kf(px) for aws process. 
The symbols represent the results from iterations for some 
cases of {po, qo, Q 2 } values studied. Different colored symbols 
represent different quark momenta. The results of variational 
method are shown by curves. The terms It. and Var. in figure 
labels refer to iterations and variational methods. 
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FIG. 6: The p± distributions Sftgi(px) (longitudinal part) of FIG. 7: The pj_ distributions 5Rgr(px) (longitudinal part) of 
bremsstrahlung photon emission. The other details are as in aws process The oter details are as in Figure. |S] 
Figure. 0] 
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FIG. 8: The peak positions of the pj_ -Kf(p±) and p± -Kgi(px) 
distributions corresponding to various processes and photon 
virtuality Q 2 and quark momenta. All data points shown are 
peak position results obtained from iterations method. The 
solid curves are the empirical power law fits shown on the fig- 
ures. One should note that the x-scale is x b T for figures (a,b) 
and scale is xi for figure (c) (for details see text). These em- 
pirical fits give a rough estimate of the optimized variational 
parameter for use in variational method. The labels b,a rep- 
resent the processes. The lables T,L represent components. 



as can be seen in Figs lbl7l However, the integrand of 
/ d 2 pj_3?<7(p^), i.e., p±dig(p±), exhibits peaking be- 
haviour. Therefore, we propose that the variational 
paramter may be fixed as the peak positions of these 
respective distributions. One should notice the x-scale is 
Xrp as defined in Eg 1451 for both transverse and longitudi- 



and aws process for Q 2 < lOGeV 2 . 
ure, the transverse distributions (i.e., 



nal distributions. The data is fitted by a function whose 
formula and the coefficients are mentioned in the figure 
(a) for both transverse and longitudinal parts. We have 
chosen a fitting function of type y = a + bx p because, as 
x — > 0, we want the peak positions to saturate as in real 
photon case (for real photons A v > 0.31). 
Figure |HJb) shows the similar results for aws process at 
Q 2 > lOGeV 2 . Here again the x scale is x\. Moreover, 
the distributions p^ • SRf(px), p±$lg(p±) have approxi- 
mately same peak position values for different po,qo, Q 2 . 
The fit functions and coefficients are shown in Figure 

EP>). 

Figure Etc) shows the peak positions for bremsstrahlung 

In this fig- 
3tf(px)) for 

bremsstrahlung and aws have approximately same peak 
positions for different po,qo, Q 2 values. Notice that the x 
scale is the x\ variable defined in Eg 1431 Surprisingly, all 
these peak position values scale with x\ variable rather 
than the usual x\ for transverse components. This might 
be indicating that for small virtuality, the relevant scale 
is xi, coinciding with real photon scale. Similarly, the 
longitudinal components p±$ig(p±) have same peak po- 
sitions for these two processes. The fit functions and co- 
efficients are shown in figure. Using the formulae given in 
Figures|5Ia,b,c) for the case of virtual photon emission, 
one may choose the variational parameter A v {po, qo,Q 2 ) 
to be around these peak position values, 
ft is interesting to compare the results for real and 
virtual photon peak positions. The peak positions of 
p± • Sif(px) distributions for real photons are given 
by empirical fits in Fig[3fc). The peak positions of 
p_i_ ■SRf(px), P±$tg(p±) distributions for virtual photons 
for Q 2 < lOGeV 2 are given in FiglHJc). In Figure E0 we 
compare these three results. As seen in figure, the peak 
positions for real photon and the transverse part of vir- 
tual photons are approximately the same. It should be 
noted that the numerical calculations have errors such 
as the iterations method has convergence errors, peak 
search has errors, errors in empirical fit of peak positions 
etc. Further, a s value used is 0.2 for real photons and 0.3 
for virtual photon studies. This close matching of peak 
positions of real and virtulal cases is rather surprising, as 
a strong Q 2 dependence is expected for virtual photons. 
Similar surprising results were already presented in pH 
regarding photon emission function. 



dp [n F (r ) - n F (p )] 
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In the previous section, we compared the results from 
variational and iterations methods and gave empirical fits 
to optimized values of variational parameters. Using the 
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FIG. 9: The peak positions of the px -Kf(px) for real photons 
represented by red circles, same as in FigEfc). The peak 
positions px-SRf(px), p±^tg(f>±) distributions corresponding 
to various processes and photon virtuality Q 2 < WGeV 2 are 
shown by black curve and green curve respectively. 



FIG. 10: QTIij plotted as a function of Q 2 /T 2 for a 
photon energy of 50GeV. The transverse components of 
bremsstrahlung, aws, the insignificant contribution from lon- 
gitudinal parts are shown. Symbols represent variational cal- 
culations. For reference the curves shown are taken from 12111. 



variational method and the variational parameters given 
in Figures|8jL,b,c, we repeated the variational calculations 
for all the values of {qo,q,Po}- These calculations cover 
transverse and longitudinal parts of bremsstrahlung and 
aws processes for which we have reference distributions 
from iterations method. Now, we compared these two 
sets of j»x distributions. It was noticed that the trans- 
verse bremsstrahlung distributions were exactly repro- 
duced for all {qo,q,po} values. The longitudinal contri- 
butions to bremsstrahlung were not not well reproduced, 
as the variational data showed oscillations around the 
iterations distributions. The transverse and longitudi- 
nal parts of aws were reasonably well reproduced for all 
Q 2 < lOOGeV" 2 . The aws distributions showed sensi- 
tivity to variational parameter. Therefore, it should be 
noted that the empirical optimized variational parame- 
ter is a good approximation, however, one needs to vary 
the variational parameter around these empirical values 
to converge the variational results. 

The imaginary part of retarded photon polarization ten- 
sor (represented by SsHr) is calculated using the p_i_ in- 
tegrated values as given in EqH3(seeEq.l6of O). The 
required px integrated values were generated using vari- 
ational method with empirical variational parameters. In 
this Eg 1461 one should note the factor 1/mjj in the lon- 
gitudinal part for reasons explained before. All terms in 



this equation contributing to QUr are calculated. Fig- 
ure El shows the QTIfl plotted as a function of Q 2 /T 2 
for a photon energy of 50GeV. The results of variational 
method are represented by symbols in the figure. The 
transverse components of bremsstrahlung and aws rep- 
resented in figure by b,T and a,T. The longitudinal con- 
tribution to aws is shown as a,L. The results from vari- 
ational method (symbols) for various components have 
been normalized to approximately match with respective 
curves at Q 2 /T 2 — 1. At this photon energy (50GeV), 
the contribution of bremsstrahlung is completely negli- 
gible. The bremsstrahlung is insignificant because the 
qo is very high. The transverse component of aws only 
contributes, with a small contribution from longitudinal 
component above Q 2 > 20GeV 2 . The curves represent 
the imaginary polarization tensor by empirical method 
proposed in |2l|. It can be seen that the variational 
has predicted reasonably well the the results o f 1211 for 
Q 2 /T 2 in the range of 1 — 250. As shown in [l!|, the 
multiple re-scatterings in the medium only marginally 
increase the ^$Hr at low Q 2 . However, the re-scattering 
corrections smooth out the discontinuity at the tree level 
threshold Q 2 = 4M^ [ig|. 
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IV. CONCLUSION 

The photon emission processes from the quark gluon 
plasma have been studied as a function of photon mass, 
considering LPM suppression effects at a fixed tempera- 
ture of the plasma. Self-consistent iterations method and 
the variational method have been used to solve the AMY 
and AGMZ integral equations. We obtained the 3?f(pj_), 
3ftflr(p_i_) distributions as a function of photon mass, pho- 
ton energy and quark momentum. The corresponding 
distributions from variational method have been com- 
pared with the results of iterations method for validating 
the variational approach. In order to fix variational pa- 
rameter, the peak positions of the pj_ distributions have 
been studied in detail for both real and virtual photons 
emission using iteration method. We identified relevant 
dynamical scales for the peak position values for real and 
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